Modifying gravity with the Aether: an alternative to Dark Matter 
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There is evidence that Newton and Einstein's theories of gravity cannot explain the dynamics 
of a universe made up solely of baryons and radiation. To be able to understand the properties of 
galaxies, clusters of galaxies and the universe on the whole it has become commonplace to invoke 
the presence of dark matter. An alternative approach is to modify the gravitational field equations 
to accommodate observations. We propose a new class of gravitational theories in which we add 
a new degree of freedom, the Aether, in the form of a vector field that is coupled covariantly, but 
non- minimally, with the space-time metric. We explore the Newtonian and non-Newtonian limits, 
discuss the conditions for these theories to be consistent and explore their effect on cosmology. 



I. INTRODUCTION 

Despite the tremendous successes of contemporary cos- 
mology, there is a nagging problem that refuses to go 
away. If we try to measure the total gravitational field in 
the universe it far surpasses what we would expect from 
the baryonic mass we can see. This is true on a wide 
range of scales. On Kiloparsec scales it is well known 
that the velocity of objects in the outer reaches of galax- 
ies are moving around the central core at much greater 
speeds than what one expect from Keplerian motion due 
to the stars and gas. On Megaparsec scales it has been es- 
tablished that the random motion of galaxies in clusters 
is too large for these systems to remain gravitationally 
bound due simply to the visible mass. And on tens to 
hundreds of Megaparsecs there is evidence for structure 
in the distribution of galaxies which should in principle 
have been erased by dissipational damping at recombi- 
nation, when the universe was a few hundred thousand 
years old. 

There is a solution to this problem. One can invoke 
the existence of an exotic form of matter that does not 
couple to light. It is cold and clumps easily to form 
bound structures. The dark matter |l[ will enhance the 
energy density of galaxies and clusters and can be mod- 
eled to fit almost all observations. It will also sustain 
gravitational potential wells through recombination and 
reinforce structure on large scales. A cosmological the- 
ory based on the existence of dark matter has emerged 
over the past twenty years with remarkable successes and 
predictive power [2|, |3j. Laboratory searches are under 
way to find tangible evidence for dark matter candidates 
which go beyond their gravitational effects. 

One can take a different point of view. At the mo- 
ment, all evidence for dark matter comes from its dy- 
namical effect on visible objects. We see dark matter 
through its gravitational field. Could it be that our un- 
derstanding of the gravitational field is lacking? This 
possibility has been mooted before. It has been pro- 
posed that the Newton-Poisson equation, V 2 $ = 4nGp 
(where $ is the gravitational potential, p is the energy 
density and G is Newton's constant) should be modi- 
fied to V • [/(|V$|/a )V$] = 4irGp where f(x) = 1 in 



the strong field regime and f(x) ~ x in the weak field 
regime. In regions of low acceleration, gravity is boosted 
above the standard Newtonian prediction and an / can 
be chosen to fit galactic rotation curves Such a the- 
ory, dubbed Modified Newtonian Dynamics (MOND) has 
proven very effective and it has recently has been pro- 
posed that such a behaviour can emerge from the low 
energy, non-relativistic limit of a fully covariant theory 
(see [5, 6] for various approaches). 

MOND is not without problems. It has been shown 
that it is less effective at resolving the missing mass prob- 
lem on the scale of clusters of galaxies. Indeed it has 
been shown by Sanders Q that to correctly account for 
the mismatch between luminous and dynamical mass in 
clusters one must invoke a small fraction of massive neu- 
trinos, with a mass of approximately 2 eV. This result 
has recently been reconfirmed with weak lensing data 
presented by Clowe et al [1] and the subsequent analysis 
by Angus et al 0. Given that neutrinos exist, are mas- 
sive and the mass required falls within the allowed range 
constrained by laboratory measurements, this solution to 
the missing mass problem in clusters is not outlandish. 

In this paper we show that it is possible to modify 
gravity by introducing a dynamical Aether (or time-like 
vector field) with non-canonical kinetic terms. Our pro- 
posal builds on the extensive analysis of Einstein- Aether 
theories undertaken by Jacobson, Mattingly, Carroll, Lim 
and collaborators [H1>E3 and follows along a long series 
of proposals by others [Ilj. As the Aether vector field has 
a non- vanishing expectation value it will dynamically se- 
lect a preferred frame at each point in spacetime (i.e. the 
frame in which the time co-ordinate basis vector dt aligns 
with the direction of the Aether field A) . This violates 
local Lorentz invariance (and gauge invariance). Conse- 
quently, Aether theories traditionally have been used as 
phenomenological probes of possible Lorentz violation in 
quantum gravity. 

As there has been recent interest in modifying gravity 
by using additional scalar and vector fields it is worth 
initially comparing our approach to other attempts. An- 
other group of theories retains local Lorentz invariance 
but introduces a vector field, the action of which breaks 
U(l) symmetry through a non-gauge invariant potential 
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or kinetic term. These theories have variously been found 
to be able to model inflation [12| and dark energy [l3j , the 
latter invoking a set of three identical vectors along mu- 
tually orthogonal spacial directions (the "cosmic triad" ) . 
Such a vector field, coupled to scalar fields, has also been 
considered in the context of alternatives to dark matter 
Quite distinctly, it has been suggested that the re- 
cent acceleration of the universe could be accounted for 
by allowing for non-canonical kinetic terms in the action 
of the electromagnetic field [HI]. Recently the possibil- 
ity of a scalar field coupled to the Aether as a cause of 
inflation has been examined (l6j . 

We will first lay out the formalism for our proposal, 
with the full field equations. We will then proceed to an- 
alyze them in the non-relativistic regime and show that 
it is possible to naturally obtain modifications to New- 
tonian gravity. The physical consistency of the theory 
is discussed in the weak field regime as are constraints 
from the Solar system. We then briefly explore the pos- 
sible impact on the expansion of the Universe, showing 
that this modification of gravity can lead to accelerated 
expansion at different stages of the evolution of the uni- 
verse. A specific proposal for such a theory is presented 
and we conclude by discussing a series of open problems. 



value for A). We allow for these different possibilities by 
deriving a general form for the field equations below. We 
will comment on these models in the discussion. 

The gravitational field equations for this theory, ob- 
tained by varying g af3 (see [TtJ but also [HI]) are 



G'T 1 i n^/^irrmiatter 
a/3 = J- a/3 + OTTLrl a £ 



i/3 



(3) 



where the stress-energy tensor for the vector field is given 
by 

f a f3 = -V a (J 7 '(J(aAp- ) -J' T ( a A l3 - ) -J( a p ) A' T )) 



J ' dIC 

-a/3 



(4) 



Brackets around indices denote symmctrization and Y a p 
is the functional derivative 

Y „ -V 4"V At 5 -^—^ 



II. THE THEORY 



The equations of motion for the vector field, obtained 
by varying A@ are 



A general action for a vector field, A coupled to gravity 
can be written in the form 



d^x. 



R 



16ttG 
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Sm 



(1) 



where g is the metric, R the Ricci scalar of that metric, 
Sm the matter action and C is constructed to by gener- 
ally covariant and local. Sm only couples to the metric, g 
and not to A. We shall use the metric signature (-,+,+,+) 
throughout. 

For most of this paper we will restrict ourselves con- 
sider a Lagrangian that only depends on covariant deriva- 
tives of A and we will consider a A that is time-like. Such 
a theory can be written in the form 



M 2 



IGttGn 
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K. = M- 2 K, a0 V a A~>V„A° 
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where Ci are dimensionless constants and M has the 
dimension of mass. A is a non-dynamical Lagrange- 
multiplier field with dimensions of mass-squared. Note 
that it is possible to construct a more complicated K, 
by including different powers in A and its derivatives. 
Indeed it is possible to show that Bekenstein's theory 
of modified gravity [j| is formally equivalent to a the- 
ory with such an extended K, (though with a more exotic 
method of achieving a non- vanishing vacuum-expectation 



V a {T'J a p )+T'y fj = 2\A P (5) 
where once again we define the functional derivative 

Variations of A will fix A^A^ = — 1. By inspection, 
transvecting both sides of ((5|) with A 13 leads to a solution 
for A in terms of the the vector field and its covariant 
derivatives. 

These equations allow us to study a general theory of 
the form presented in equation [T] with a time-like vector 
field. For our particular, restricted choice of K, we have 
Y a p = -a [(VA)(VM ? ) - {V a A v ){VpA v )) and y p = 
0. 



III. THE NON-RELATIVISTIC REGIME: 
NEWTONIAN AND MONDIAN LIMIT 

Having established our general theory, we can now ex- 
plore its properties in various different regimes. We start 
in the static, weak field, non-relativistic limit. We must 
expand both the metric and vector field around a fixed, 
Minkowski space background. The constraint equation 
very clearly fixes the zeroth order contribution of the 



vector field to be A^ 

9af3 

A a 



Sq. Hence we have 
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We will only consider terms which are linear in our pa- 
rameter e and restrict ourselves to the Poisson gauge (so 



too 



-2$,hij 



-2^6a), We have that 



VA = eS + 0(e 2 ) 



(G) 



where = 0, S 0m = ^$ and S i0 = + B°) and 
(i) indexes the three coordinates, (x 1 , x 2 , a; 3 ). Terms in 
VAVA and B may be consistently regarded as being of 
0(e 2 ) or greater and so can be ignored in what follows. 
An immediate consequence of this is that, to first order, 
the ijth component of the Aether's stress energy tensor 
disappears. The ijth component of the Einstein equa- 
tions then yields $ = ^. The equality of scalar potentials 
is true of General Relativity in the absence of anisotropic 
stress and simplifies the equations dramatically. 

The modified Newton-Poisson and vector equations 
then are 

2V 2 $- (ci - c 3 )V .{T'V B°) - A = 8nGp 
2c 3 V.(J c ''V$) + 2c 1 V.(J c ' , V($ + S )) = -2A (7) 



The constraint fixes K, to be 



K, = — ci 



(V$) 2 
M 2 



(8) 



Taking c\ < ensures that K, is positive. We can manip- 
ulate the system of equations (O so as to obtain the field 
equation for <i>: 



V.((2 + ciJ r ')V$) = 8nGp 



(9) 



If such a theory is to have V.(|V$| V$) cx p in the 
limit of small V$ we require that: 



lim \2 + c\T'\ cx /C' 

|V$|<<A/ 



(10) 



Integrating then we have that J 7 = a/C + /3/C' where 
a and /3 are constants. Hence we can construct a theory 
with 'MONDian' limit on galactic scales, this holding as 
long as we identify M with something of the order of ao 

The limit as K, — > is worth considering in more de- 
tail. Considering a test particle a distance r away from an 
isolated source, the Modified Poisson Equation dictates 
that for |V$| << M the gravitational force will vary as 
1/r. This combined with the geodesic equation suggests 
that in this regime the metric component goo grows ap- 
proximately as + ln(r) and therefore becomes ill defined 
in the limit of r — > oo. This has been a generic feature 
of past metric theories of MOND. It was noted in [26[ 
that the value of this growing term will typically vary by 
around only 10~ 6 from the radius of onset of MOND in a 
system to the present Hubble radius (> 10 27 cm). There- 
fore, the weak field approximation is unlikely to be seen 
to break down. The situation becomes more complicated 
if the theory is such that MOND only arises in a cosmo- 
logical background (for instance if A 2 is not fixed and 



the 'scale' M is contingent on its value at the particular 
cosmological era) . It would only be consistent then to for- 
mulate the weak field limit in a FRW background rather 
than Minkowski spacetime and the asymptotic form of 
the metric may be expected to differ. 

The Solar System can supply us with stringent con- 
straints on the weak field limit of these theories. Accel- 
erations are typically substantially larger than M (again 
assumed to be ~ ao) and so requiring concordance with 
observations could constrain the possible form of T for 
large JC. Poisson's equation V 2 $ = AnGp is an excellent 
approximation in the Solar System and thus we shall ex- 
pect the contribution of c\T' in @ to be small. We may 
expect then that T' in this limit can be expanded as a 
power series in inverse powers of /C 5 . That is: 



lim ,ra=V^ 



|V*|>>M 



AC 3 



(11) 



where r\i are constants. 

Consider the leading term 771//CX We would expect 
such a term in spherical symmetry to result in a constant 

The existence 



anomalous acceleration equal to V c ^i M _ ^ 



of such a term is particularly constrained by observed 
bounds on the variation of Kepler's constant GMq. For 
instance, observations between Earth and Mars restrict 
the acceleration to be less than approximately 10~ 9 ms -1 
[20l | . We shall see that 77 and c values are typically of or- 
der unity so the above is rather restrictive. Even under 
the assumption of spherical symmetry, the field equations 
for the theory are enormously complicated. The inher- 
ent nonlinearity in the weak field limit provided by T' 
presents a considerable challenge 21]. We note however 
that the fixed norm constraint on A will, in the weak field 
limit, force terms of the form VAVA and VVA to be at 
most of the order of terms comprising the components of 
the Einstein tensor G. For the limiting form of T 1 given 
in (jlip we expect terms in the Aether stress energy tensor 
to be schematically of the form -^cj [V V AandV AVA] 

and higher order derivative terms suppressed relatively 
by factors of . 

At Mercury, the ratio |V$|/M is of order 10 s . Provi- 
sionally neglecting terms in T' of K.~ 3 (see above), we 
then expect corrections to terms in G in the inner solar 
system to be of order 10~ 16 CjCti/ck- It is tempting to con- 
clude that two 'Parameterized Post Newtonian' (PPN) 
parameters measurable by inner solar system effects, (3 
and 7, would generally be expected to deviate from the 
predictions of General Relativity by a similar order. The 
complete set of PPN coefficients have been obtained for 
the case T{JC) cx K, ([Iol|. 3rd reference) . In particular 
the coefficients describing 'preferred frame' effects were 
found to be only consistent with experimental bounds for 
specific combinations of the c^; they are expected to be 
a particularly strong test of more general forms of T . 

Additionally we note that the asymptotic behaviour of 
T is consistent with our assumption that the term M 2 ^-g 
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is second order in perturbations. We see that for JC << 1, 
M 2 T -> |V$| 2 , and for JC » 1, M 2 T shall be at most 
aiM|V$|. 

Finally we briefly consider the effect of allowing for 
first-order spatial components of A, A 1 . It has been 
found [23] that in Bekenstein's theory of modified gravity 
the growth of large scale structure is necessarily accom- 
panied by growth in A 1 so it is not unreasonable to expect 
that such behaviour shall be present in the model con- 
sidered here. It may be readily checked that in the static 
weak-field limit that an A 1 of order e will only contribute 
to Too and Ta at order two in e and above. However, as 
the numerator in JC is a second order quantity, it shall 
generally be affected. If A % has only a radial component 
(i.e. A r ^ 0), it may be checked that JC is modified as: 



Cherenkov radiation. This radiation is expected when 
gravitational waves propagate subluminally. With the 
Aether, the usual transverse traceless modes exist along 
with three coupled Aether-metric modes. The speeds of 
propagation of each mode are functions of the c$ and have 
been calculated in [Toj for perturbations in Minkowski 
spacetime. The squared speed s 2 t of the transverse- 
traceless mode is: 

4 = i/(i + ( Cl + C3 )) 

and the squared speeds of the transverse Aether and 
trace Aether-metric modes (sf A and s\ r respectively) are 
given by: 



jc = 



-ci(V$) 2 + c 2 (V r A'-) s 
M 2 



(12) 



where the covariant derivative V ' r A r is a of order two 
or greater in e. We shall see later that c 2 is preferably 
the same sign as c\ . The effect of gradients of the radial 
component of the vector field then is to decrease JC for a 
given |V$|. Recall that the onset of MONDian behaviour 
coincides with JC << 1. Therefore an A r (r) in this model 
will generally hasten the onset of this limit, the effect 
being to further increase the gravitational field for a given 
p (see ©)• It requires further work to see whether this 
can appreciably counter MOND's problems on the scale 
of clusters of galaxies. 



IV. FURTHER CONSTRAINTS 

Recall that T tends to alC + (3JCi for small JC (i.e. 
far from a source). Therefore, in general, when consider- 
ing classical perturbations one is effectively considering a 
theory with J- ~ JC, as in the Einstein- Aether theory with 
minimal couplings. This can be used to study the consis- 
tency of these theories in the perturbative regime. Lim 
[24| has considered the dominant term in the limit where 
metric and vector field perturbations decouple. The vec- 
tor field propagates in flat spacetime and allows for a de- 
composition of perturbations into spin-0 and spin-1 com- 
ponents. The requirement that the perturbations can be 
consistently quantized, that the spin-0 component propa- 
gates subluminally and nontachyonically when quantized, 
and that metric perturbations do not propagate superlu- 
minally place the following restrictions on the Cf. 



Mi 



-cia)) 
(<Wci)(2- 



[2(1-C 2 ) 2 +Ci23(l-C 2 -Ci23)] 



where 



■ijk. 



— c l + cu + c k 



and Ci are the opposite 



sign to those considered in [lj| ■ 

Computing the expected degree of gravitational 
Cherenkov radiation for the above modes and compar- 
ing this with observational bounds, it has been found 
25| that for the case !F(JC) = JC that extremely stringent 
constraints could be placed on the Ci. For instance, they 
may be satisfied when the Q are mutually related such 
that the metric and Aether-metric modes propagate at 
precisely the speed of light. If this is not the case then 
the magnitudes of the Cj must be severely diminished. It 
is expected that such an analysis could powerfully con- 
strain more general combinations of T,c%, though any 
nonlinearity is expected to complicate the results when 
propagation is considered against a curved background 
spacetime. 

Note that such an analysis is by no means complete. It 
has been shown that a very restricted class of Einstein- 
Aether theories do not have positive Hamiltonians and 
therefore are inherently unstable at both the classical and 
quantum levels [27| . Furthermore we are considering non- 
linear functions of JC and hence instabilities may arise 
in non trivial backgrounds. A more detailed analysis of 
individual cases for T is needed yet first indications are 
that these theories are healthy. 



COSMOLOGY 



ci < (13) 
c 2 < (14) 
ci + c 2 + c 3 < 0. (15) 

Additional constraints can be obtained via astropar- 
ticle physics. The observation of ultra-high energy cos- 
mic rays implies a lack of energy loss via gravitational 



This class of modified theories are generally covariant. 
This gives us the possibility of exploring their properties 
on large scales and in particular we can consider the case 
of a homogeneous and isotropic universe in which the 
metric is of the form ds 2 = —dt 2 + a(t) 2 dx 2 where t is 
physical time and a(t) is the scale factor. The vector 
field must respect the spatial homogeneity and isotropy 
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of the system and so will only have a non-vanishing T 
component; the constraint fixes A = (1,0,0,0). The 
energy-momentum tensor of the matter is of the form 
rpmatter = p U a U () + P{g a p + U a U p) where p is the energy 
density P is pressure and we have introduce a four- vector 
U which satisfies g a bU a U b = — 1. 

The equations of motion simplify dramatically with 
these symmetries and we find that: 



JC 



3H 
aH 2 



M 2 



does introduce an undesirable discontinuity in T at the 
origin; however, as the behaviour in (jl8p will lead to a 
constant JC (as H = 0),this form of T need only hold 
in the domain |/C/3/3| >~ Hq and does not inform the 
desirable choice of T around JC = 0. 

Accelerate expansion can be reached in a more natural 
way by positing 

for a certain range of JC For arbitrary n Equation |T 
becomes 



where H = —, the dot denotes differentiation with re- 
spect to t, and, following [l7j . we define a — c± + 3c2 + c^. 

Note that now JC is negative, unlike the nonrelativistic 
limit encountered above. This means that the dynam- 
ics of static, spherically symmetric systems on the one 
hand and of relativistic cosmologies on the other probe 
completely different branches of T. 

The modified Einstein's equations now become: 



1 + e 



2(n-l) 



H z 



8ttG 



where e = — (1 — 2n)j(— 3a) n /6.For an appropriate choice 
of 7 and n we have e < and we find that H tends to an 
attractor 



M(- 



> 1/2(1- 



[1 - T'a}H 2 



-TM = 



8nG This regime is approached asymptotically with 



-P 



[1 - 2T'a]H 2 - 2[1 - -J 7 ' a]- + T'aH - -TM 2 = 8nGP 
2 a 2 



H — H. 



<■</ 



H, 



eq 



8nGp 



These can in fact be rewritten in a more useful form: i.e. as a 3 in the matter era. This gives us a particularly 



1 - aJC 1 ' 2 



d 



T 



die Kjc 1 / 2 



H 2 



8nG 



(16) 



— (-2H + T'aH) = 8ttG( P + P) (17) 
at 

Once again these equations are general but we can see 
in their structure, interesting possibilities. If we take T = 
0, we recover the standard cosmology. We can do this by 
either setting it to or choosing JC to have c\ = — c 3 
and C2 = so that it becomes the Maxwell tensor. We 
then have a theory which will modify gravity on galactic 
and super-galactic scales but leaves the expansion of the 
universe unchanged. 

It is interesting to explore the possibility that the vec- 
tor field may affect the late time expansion of the Uni- 
verse. From Einstein's equations, we can see that for it 
to behave like a cosmological constant, A we must have 



(18) 



elegant, dynamical mechanism for approaching acceler- 
ated expansion without invoking a cosmological constant. 
A small scale is still invoked but it is naturally related to 
the acceleration scale needed to trigger the onset of the 
MONDian regime on galactic scales. 



VI. DISCUSSION 

The discussion throughout this letter has been general. 
We have neither chosen a specific form of JC or T although 
we have constrained the asymptotic form of the latter. 
It is instructive to pick a simple example. If we choose 
ci = —1 ,C2 =0 and C3 = 1 we will recover the canonical 
form JC oc F a pF a P where F is the field-strength of the 
four vector A. A possible functional form for T is then: 

Tip, = 4(/d - ln(/C5 + 1)) 

From which we recover the modified Newton-Poisson 
equation: 



which we can solve to find: T = ^(-/C) 1 / 2 - Jjj-. 

Where f3 is a constant of integration. Now note that for 
positive arguments of T we have been equating the mass 
scale M with the acceleration scale ao ~ 10~ 10 ms~ 2 ~ 
cHq where c is the speed of light and Hq is the Hubble 
constant. If we assume that is also true for negative 
values, we find that T is of order unity. I.e. there is 
a natural relation between the scale of the cosmological 
constant and the fundamental mass scale in our theory. 
A choice of such asymptotic behaviour close to JC = 



V.(- 



]V*| 

M 



M 



-V$) = AirGp 



1 



(19) 



With an appropriate choice of the value of M we recover 
the field equation considered in [22j , which was found to 
give a satisfactory fit to the terminal velocity curve of 
the Milky Way. As noted, the choice ci = — C3, C2 = as 
here forces JC to vanish in the case of spatial homogeneity 
and isotropy. Therefore with this choice of T the Aether 
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has no influence on cosmological background evolution. 
The inclusion of a nonzero C2 modifies the kinetic term 
of the metric and may have significant effects both in the 
static, weak-field limit and the cosmological background. 

There arc a number of theoretical and phcnomcnolog- 
ical issues that remain to be addressed. 

For a start, the presence of the non-dynamical 
Lagrange-multiplier field A in the Lagrangian @ is per- 
haps unappealing. Its sole role is to impose the constraint 
that is unit-timelike. The same can be accomplished 
by replacing the A-term in the Lagrangian by a potential 
term — V(A fl ) ) by dint of which at low energy acquires 
a vacuum expectation value such that A^A^ = — 1. For 
example V(A^) = i/i 4 (A^A^ + l) 2 , with fi a constant 
with dimensions of mass. At energy scales below //, one 
expects that indeed A 2 = —1, and that "radial" excita- 
tion of A (i.e. of a = fj,A) will have positive mass-squared 
m 2 . ~ /i 2 . So long as fj, 2 is sufficiently large, the low- 
energy phenomenology of this model should be identical 
to that of |2]). A more exotic possibility is to construct 
more complex K, which have minima for time-like A. In- 
deed the theory proposed in || is of this form [28[ . 

Closely related to the previous point is the fundamen- 
tal origin of such an Aether. We have kept the discussion 
general in the hope that a more fundamental theory of 
fields or strings may pin down the form of K, and T . In- 
deed, the possibility of Lorentz violating vector fields has 
cropped up in attempts to extend the standard model of 
particle physics. Most notably it has been argued that 



such a field may arise in higher dimensional theories as a 
low energy by product of string theory [301 ] . It would be 
interesting to explore the range of current candidates for 
low energy string theory to find a possible candidate for 
such an Aether field. What is clear is that a Lagrangian 
of the form we require will not appear from the standard 
perturbative approach to constructing effective field the- 
ories as in [191 ]. Non-perturbative effects must come into 
play. 

We have endeavoured to explore some of the observa- 
tional constraints. Preliminary indications are that these 
theories are compatible with Solar System constraints. 
A more detailed analysis is needed with complete spher- 
ically symmetric solutions that can then be used to cal- 
culate the Post Newtonian Parameters [3l| . We have 
also shown that, as yet there is sufficient freedom in the 
choice of T to obtain different cosmological behaviours, 
from no effect to early or late time acceleration. The next 
step is to follow in the footsteps of [HJ and calculate the 
evolution of linear perturbations. A priori it is unclear 
whether perturbations in the Aether will have the same 
effect as dark matter in sustaining perturbations through 
the Silk damping regime during recombination [33l |. In- 
deed this may be the most stringent test such theories 
have to pass to be viable alternatives to dark matter. 
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